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Abstract 

We study a semiclassical inverse spectral problem based on the spectral asymp- 
totics of [llj . which apply to small non-selfadjoint perturbations of selfadjoint h- 
pseudodifferential operators in dimension 2. The eigenvalues in a suitable complex 
window have an expansion in terms of a quantum Birkhoff normal form for the opera- 
tor near several Lagrangian tori which are invariant under the classical dynamics and 
satisfy a Diophantine condition. In this work we prove that the normal form near a 
single Diophantine torus is uniquely determined by the associated eigenvalues, possibly 
up to some natural types of symmetry which are explained in the main result. We also 
discuss the normalization procedure and symmetries of the quantum Birkhoff normal 
form near a Diophantine torus. 



1 Introduction 



Let M denote either or a compact, real analytic manifold of dimension 2, and let M 
denote a complexification of M, which is in the Euclidean case, and a Grauert tube of 
M in the compact analytic 2-manifold case. 

We study operators of the form = P + ieQ, where P, Q are analytic /i-pseudodifferential 
operators on M with principal symbols p, q, respectively, and P is self adjoint. The principal 
symbol of Ps is then p + ieq, where p is real. We will also make a non- degeneracy assumption 
on Reg, but we do not require Q to be selfadjoint. 

♦Partially supported by NSF grant EMSW21-RTG. 
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Consider a Lagrangian torus A, contained in an energy surface p~^{0) C T*M, which is 
invariant with respect to the Hamilton flow of p and satisfies a Diophantine condition (see 
section [2^ . For simplicity, we will assume: 

The Hamilton flow of p is completely integrable in a neighborhood of p~'^{0)- (1) 

This implies that in a neighborhood of A the energy surface is foliated by invariant La- 
grangian tori. 

In action-angle coordinates {x,C,) such that A is the set {C, = 0}, let 

{m) = j^2 j^^^^^^Odx (2) 

denote the spatial average of (g) . Here we take x to be a multi- valued function whose gradient 
is 27rZ^-periodic. We assume that dp and dRe{q) are linearly independent along A. Then, as 
is explained in section |3l we may make a sequence of changes of variables which transforms 
the symbol of into a Birkhoff normal form, which in this context means an asymptotic 
expansion in (^, e, h) that is independent of x to high order. Formally we may carry out this 
procedure on the level of operators by conjugating by a sequence of appropriately defined 
Fourier integral operator, obtaining what we call a quantum Birkhoff normal form for P^. 

Under some further technical assumptions which we will explain later on, one of the main 
theorems of [H] establishes, for any 5 > 0, asymptotics for the eigenvalues of P^ in an 
/i-dependent "window" 

D = {zeC : \Rez\<—, \ lmz - eReF\ < —}. (3) 

Here F = {q)\A = {q){0), and we assume that this average is not shared by any other 
invariant torus. The expansions are given in terms of a Bohr-Sommerfeld type condition and 
the Birkhoff normal form of Pg near A. 

Our goal in this work is to solve the semiclassical inverse problem of determining the 
quantum Birkhoff normal form of P(. from the eigenvalues in D, assuming the unperturbed 
operator P is known. 

Inverse spectral problems have been studied for many years, as surveyed for example by 
Zelditch |T6]. Recently, semiclassical inverse spectral problems have been investigated by 
several authors, such as Colin de Verdiere [2] [3], Colin de Verdiere-Guillemin |1] [5], Guillemin- 
Paul-Uribe |8], Guillemin-Paul [7], Guillemin-Uribe [9], lantchenko-Sjostrand-Zworski [T3] . 
Via Ngoc [15]. Often in inverse spectral problems one studies the wave trace, in the spirit 
of Guillemin [6]. The dimension 2 case is special because the eigenvalues have an explicit 
description in terms of the Birkhoff normal form and Bohr-Sommerfeld type rules, so here 
it seems most natural to recover the normal form directly from eigenvalue asymptotics. Our 
approach is taken very much in the spirit of Colin de Verdiere [3]. 

According to [H], the eigenvalues in the window (j3]) form a distorted lattice, with horizontal 
spacing ~ h and vertical spacing ~ eh. The window is of size by eh^, for some < S « 1, 
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which means that the asymptotic expansions are vahd for a comparatively large number of 
eigenvalues (on the order of as — 0. In addition, in the semiclassical setting, 

we assume we know the eigenvalues for each sufficiently small value of the semiclassical 
parameter (or for a sequence of values of h tending to 0). This provides a rich data set, from 
which we will recover information about the Birkhoff normal form using elementary order of 
magnitude arguments. 

In the perturbative setting there is also the paremeter e to consider, since the results of 
[TTj apply to all values of e in a range [h^,h^]. We will treat the cases where we are and are 
not allowed to vary e independently of h. 

Our main result, stated informally, is the following: 

Theorem 1. The eigenvalues of in D determine the quantum Birkhoff normal form of 
Pi; near A, up to some ambiguity when e is comparable to a rational power of h. 

This is stated more precisely in Theorem [121 

The plan of the paper is as follows. In section [2l we recall the setting and technical 
assumptions of [11] needed to apply one of the main results of that paper. In section [3l 
we review the normal form construction in the present context, using some notation and 
methods due to S. Vii Ngoc, and discuss symmetries and uniqueness of the normal form. In 
section HI we recall the spectral asymptotics of [H] in a precise form, which is the basis of 
the inverse result. Finally, in section [5], we prove our main result. 



2 Assumptions 

We now recall the setting of [TT]. We will state our assumptions along the lines of section 
7 of that work, in particular restricting our attention to the completely integrable case ([1]) 
rather than the most general case considered elsewhere. 

In fact, we can imagine that our main result holds (with exactly the same proof) more 
generally. In principle, one only needs that the asymptotic expansions given in Theorem 1.1 
of [11] are valid, as well as the Diophantine assumption, which is explained below in section 

2.1 Analyticity and general assumptions 

Let us assume that Pe = P + ieQ satisfies the same general assumptions as operators studied 
in [TT], which we recall here for convenience. 

When M = R^, assume that Pe = Ps{x, hD^, h) = P{x, HD^, h) + ieQ{x, HD^, h) is the 
Weyl quantization of some total symbol, which we also denote P^ = P + leQ. Assume that 
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is a smooth mapping of (e, h) into the space of holomorphic functions of (x, ^) in a small 
complex neighborhood of R^. Assume that 

\P,{x,^,h)\<O{l)g{Re{x,0) 
in this neighborhood, where g is an order function in the sense that 

g{X) <C{X- Yf' g{Y), OO, No>0, X,Y e R^ 
and that p satisfies an ellipticity condition 

When M is a compact, real analytic 2- manifold, assume that in any choice of local coordi- 
nates Pe is a differential operator with analytic coefficients, which themselves have asymptotic 
expansions in powers of h and are smooth in e. Assume also that p satisfies an ellipticity 
condition near infinity, 

b(x,oi>^(er, (x,oeT*M, iei>c, 

where we implicitly assume that M has been equipped with a Riemannian metric, so that 
the quantity (^) = ^^1 + makes sense. Assume also that the underlying Hilbert space is 
L^(M, where /i is some positive real-analytic density on M. 

The above assumptions imply that has discrete spectrum in a fixed neighborhood of 
G C for sufficiently small h, e. 

2.2 Assumptions on the classical dynamics 

Assume the energy surface p~^(0) fl T*M is non-critical, i.e. dp ^ Q along this set. For 
simplicity, we assume p~^(0) is connected. Let 

dp d dp d 

denote the Hamilton vector field of p (in local coordinates). 

By the complete integrability assumption ([T]), there exists an analytic, real- valued function 
p such that = {p, p} = with dp and dp linearly independent almost everywhere. Here, 
{■, ■} denotes the Poisson bracket. Then the energy surface decomposes as a disjoint union 
of compact, connected ifp-invariant sets, which we assume has the structure of a graph, in 
which edges correspond to families of invariant Lagrangian tori and vertices correspond to 
singular invariant sets. 

Near an invariant torus A we have real analytic action-angle coordinates {x,C,) such that 
Hp = a ■ dx for some frequency vector a when restricted to A. We refer to A as a rational. 
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irrational, or Diophantine torus if the vector a has the corresponding property. Below, we 
will consider a Diophantine torus, i.e. one such that the frequencies a satisfy 

for some Co, Nq. (Here, n = 2.) 
In action-angle coordinates near any such A, the principal symbol p takes the form 

p=p{o = ^-^+o{e)- (5) 

In particular, it is independent of x, and thus already in Birkhoff normal form. 

Recall that (q) |a denotes the average ([2]) of g with respect to the natural smooth measure 
along A, and (q) in general denotes the average in of q in the x-direction. We assume that 
(q) is not constant on any of the aforementioned "edges" consisting of families of invariant 
tori. 

Let {q)rp denote the symmetric time T average of q along the Hp-How. 

1 /■'^/^ 

(9)t (3^' = 7^ / qoexp{sHp){x,^)ds. 

^ J-T/2 

For each invariant torus A, define the interval 



goo(A) 



lim infRe(g)7^, lim supRe(g)r 

T-j-oo A T-s>oo A 



(6) 



As in [11], we have that specPg fl {| Re 2;| < 5} is contained in a band 

Im^; 



inf Qoo(A) - o(l), sup goo(A) + o(l) 

A A 



as e, /i, 5 — 7- 0. 



From now on, fix a single Diophantine invariant Lagrangian torus A, set F = (q) 
assume that dp and dlie{q) are linearly independent along A. 



and 



With all the assumptions above, and in particular assuming complete integrability, the 
last global assumption needed is that the averages satisfy 

ReF^g^(A'), AVA- (7) 
Without assuming complete integrability, a different assumption is needed (see [H], (1.24)). 

The eigenvalue asymptotics of [11] are valid in the {h, £)-dependent rectangle ([3]) for suffi- 
ciently small h and assuming e G [h'^'^,h^], where M is a fixed integer, which can be chosen 
arbitrarily large, and 5 > is also fixed but can be chosen to be arbitrarily small. 

Remark 2. The results of [H] apply also to a finite collection of Diophantine tori sharing the 
value F, in which case the set of eigenvalues in ([3]) is simply the union of the contributions 
from each individual torus. We will not consider the problem of separating the contributions 
of several tori. The global condition ([7]) implies the value F is unique to A. 
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3 Quantum BirkhofF normal form 



In this section we present the quantum Birkhoff normal form construction near a Diophantine 
torus for a perturbed symbol, and discuss issues of uniqueness for the normal form and 
normalizing change of variables. The presentation in section 13.11 follows an unpublished 
note by S. Vu Ngoc, and we are grateful to him for allowing us to reproduce it here. For 
a similar exposition in the case of a local minimum, one could read [1]. Though we only 
need to consider dimension 2, it is natural to carry out the discussion in dimension ra, as no 
changes are needed. We will work on T*T"', assuming we are in a microlocal model where 
the Diophantine torus in question corresponds to the section = 0}. 

3.1 Normal form construction 

Let us identify symbols on T*T" with their formal Weyl quantizations. The Moyal formula 
defines a product operation on symbols which corresponds to composition of the correspond- 
ing operators. We denote by [■, ■] the associated bracket operation, which on the level of Weyl 
quantizations is simply the commutator bracket. 

The normal form construction may be summarized as follows: We make a sequence of 
analytic, symplectic changes of variables, which transform to a symbol which is indepen- 
dent of X to higher and higher order in h). On the level of operators this is formally 
equivalent to conjugating by a sequence of Fourier integral operators. The resulting sequence 
of symbols is convergent in the space of formal power series in e, h). We define the quan- 
tum Birkhoff normal form (QBNF) of near the Diophantine torus A to be the formal 
limit of this procedure, while if we truncate the procedure after finitely many steps we get 
a well-defined analytic change of variables. As we remark in section 13.21 the normal form is 
well-defined once the Diophantine frequencies a are fixed. 

Later, we will often write the QBNF as a formal expansion 

where Pjki is a homogeneous polynomial of degree j in ^, and p(°°) denotes the entire formal 
expansion. 

For the moment, however, it is convenient to use different notation. Consider a grading 
in {^,e, h) which counts the power in C, plus twice the power in [e, h). Let 0{N) denote the 
associated order classes. Here we do not attach any special significance to the number two 
(e.g. three would work as well), but we note the convenience of this sort of grading: Because 
the Moyal formula has an asymptotic expansion in powers of (f^, the higher weight of 
h ensures that each time we lose an order in ^ we gain one in h. This implies that the main 
contribution in the bracket ^f-, ■] comes from the Poisson bracket of the two symbols. 
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When Kj = 0{j) and Ki = 0{i), their Poisson bracket satisfies 



k 



By what we have said above, we see that ji[Kj, K(\ = {Kj, K^} + 0{j + i). 

Proposition 3. Suppose that P = Pi + 0{2) is analytic in x and where Pi = a ■ C,, and 
that a satisfies the Diophantine condition Q). Then for all N > 1 there exist functions 

G'(^) = Gi + + ■ ■ ■ + Gjv, P^""^ = Pi + P2 + --- + Pn, Rn+i, 

which are analytic in x, with Gj, Pj and Rj homogeneous of degree j with respect to the 
grading described above (thus polynomials in C,), such that 

exp (f adG(iv) ) P = P(^) + Rn+i + 0{N + 2) (8) 

with each Pj independent of x. 

Here, adcP = and exp(f adG)P = exp(f G)Pexp(-;iG). 

Remark 4. Note that although it is not ruled out by the notation, it will follow from the 
proof that no half-powers of /i or e appear in the normal form. 

Proof. We proceed by induction on the order A^. By assumption the claim holds for = 1, 
with Gi = 1, and R2 representing the homogeneous terms of degree 2. 

Inductively if (ED holds, then setting G^^+i^ = G^^) + G;v+i, for some function G^v+i — 
0{N + 1) to be determined, we claim that the only new term modulo 0{N + 2) is given by 
{Gn+i, Pi}- Indeed, by the Campbell-Hausdorff formula 

exp adciiv+i)) P = exp (| ^^Gm+i +i adc-civ)) P 

= exp [l ads) ° exp adc^+J o exp (f adgciv)) P 

where S = C(|[Gw+i, G^^)]) = 0{{N + 1) + 1 - 1) = 0{N + 1). Meanwhile 

exp adc^+i) o exp (f adG(iv)) P = exp adc^+J (P^^^ + Rn+i) 

= 0{N). 

Thus applying exp(^ ads) only affects the terms of order 0{N + (A^ + 1) — 1) = 0{2N) = 
0{N + 2) (as N >2). Therefore we have that 

exp (i adG(iv+i)) P = exp adc^+J (P(^) + P^+i) + 0{N + 2) 

= P(^) + P;v+i + i adG,^,(P(^) + P^+i) + 0{N + 2) 
= P(^) + Pjv+i + i adG,^,(Pi + 0{2)) + 0{N + 2) 
= P(^) + P;v+i + i adc^^, Pi + 0{N + 2), 
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Because the bracket [■, ■] reduces to the Poisson bracket when one of the arguments is at 
most quadratic, we have ^ ad^jv+i Pi = jiIGn+i, Pi] = {Gn+i^ Pi}, so 

exp adG(iv+i)) P = + Rn+1 + {Gtv+i, Pi} + 0{N + 2), 

as claimed. 

To make the homogeneous order + 1 terms independent of x, it suffices to solve the 
"cohomological equation," 

{Gn+1, Pi} = (Rn+i) — Rn+1, 

for Gn+1, where (Rn+i) is the x-average of -Rat+i as in ([21) (thus -Rat+i and {Rn+i) are 
cohomologous) . We then set Pn+i = (Rn+i) and let -RAr+2 represent the homogenous order 
+ 2 part of the 0{N + 2) error terms, which are analytic. 

Because Pi = a ■ ^, we have {Gn+i, Pi} = —Hp^Gn+i = — (a ■ dx)GN+i- Expanding Gn+i 
and Rn+1 in Fourier series 

GN+i{x,i)= ^ G;v+i(n,Oe^"" 
Rn+i{x,0= 5^P^+i(n,0e^'^^ 

we have 

(a ■ dx)GN+i = ^{ia- n)G^+i(n, Oe'""". 

neZ2 

When n E Z"', because a ■ n does not vanish by (jl]), we may set GN+i{n,^) = 
—iRN+i{n,^)/{a ■ n), to obtain (a ■ dx)GN+i = Rn+i — Rn+i{0,O = Rn+i — (Rn+i) and 
thus solve the cohomological equation. Furthermore, by the Diophantine condition (jlj), we 
have 

\GN+i{n,0\ = < |nro|P;v+i(n,OI, 

t\a ■ n\ 

and because Rn+i is analytic in x in a (complex) neighborhood of T" x {0} C (C"/27rZ") x 
C", so too is Gn+1- □ 

Remark 5. As mentioned above, exp(|- adc.(iv)) formally represents conjugation by a mi- 
crolocally defined Fourier integral operator exp(^G'^^^), and such conjugation microlocally 
implements the symplectic transformation expiJg(jv). To define such operators globally, by 
means of a suitable weight function, requires care because the Hamilton fiow is complex, and 
one must work on complex deformations of T*M. The details are explained in jllj . 
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3.2 Symmetries and uniqueness of the normal form 



In this section we discuss symmetries and uniqueness of the Birkhoff normal form. 

We first remark that there is always some flexibility in choosing action-angle variables 
(x, ^) near an invariant torus (here x G R"'/27rZ" represents the angle variables, and ^ G R"" 
the action variables). If A G GL(n, Z) and ip is any smooth function on R", then 

K : {y, T]) ^ (x, = (A^'y + ^^^(r/), A't]) (9) 

gives a well-defined smooth, symplectic change of variables (which is analytic if ijj is analytic) 
on R"/27rZ" x R", and thus a new set of action-angle coordinates {y,f])- 

This transformation also preserves independence of the angle coordinate, and thus takes 
one asymptotic expansion which is in normal form (to order A^) to another. More precisely, 
if 

p{x, e, e, h) = p(e, £, h) + 0{N) = a ■ e + 0(2), 
then in the new coordinates 

{p o K){y, r], e, h) = p{A'ri, e, h) + 0{N) = (Aa) ■ 7] + C(2). 

Once we fix a choice of frequencies a, though, A must be the identity (because of the 
Diophantine assumption) , which means we only have maps of the form 

{y,r])^{y + Vi>,r]), (10) 

which do not affect a normal form expansion because they are independent of x. 

Our aim is to show that the formula fIlOp gives all transformations which preserve indepen- 
dence of the angle variables in a function or asymptotic expansion, while not affecting the 
frequencies a. Appendix A.l of Hofer-Zehnder |T2] essentially contains a proof using gener- 
ating functions that if a real symplectic diffeomorphism [y, rf) i— )■ (x, ^) satisfies ^ = h{rj) with 
det(fe^) 7^ 0, then the mapping is of the form ([9]). The only difference is that the argument 
given there is local, and they end up with a more general type of transformation. In our 
case, it turns out one can make the formula apply globally, and then because we are on a 
torus periodicity forces the simpler form ([9]), which reduces to (fTO!) assuming a is unchanged. 
This argument will be given in Proposition [7] below. 

Before proceeding, however, we note that because we have used complex symplectic trans- 
formations in our reduction to the normal form, it is natural to also consider symplectic 
biholomorphisms in a small complex neighborhood of T" x {0}, for example allowing ip to 
be complex- valued in (fTOj) . Here when we say a transformation is symplectic or canonical we 
mean that the mapping preserves the standard symplectic form a, which is a form of type 
(2, 0), given in coordinates by 

n 

a = dC^j A dxj. 
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The standard fact that any symplectomorphism locally admits a generating function carries 
over to the complex setting: 

Proposition 6. If n : {y,f]) h-> (x, ^) = {a{y,ri),b{y,r])) is biholomorphism between neigh- 
borhoods U,V ^ C^" and det{hri) 7^ 0, then there exists a holomorphic function (fiiy,^) such 
that 

K : {y,dyip) ^ (<9c<^,0- 

The proof is exactly the same as in the real case (see appendix A.l of Hofer-Zehnder [12] 
for example), if we note that the implicit function theorem holds for complex maps, and 
where normally we use Poincare's lemma we instead use the Dolbeault-Grothendieck lemma 
(for d instead of d, as normally stated). 

Using this fact, we may now prove 

Proposition 7. Suppose k: {y,ri) 1— )■ (x,^) is a symplectic biholomorphism between open 
sets U,V C (C"/27rZ") x C", where U and V are small neighborhoods of (R"/27rZ") x {0}. 
Suppose C, = b{ri), where b is a biholomorphism defined near G C", such that b{0) = 0, 
bjj{0) = 1. Then in a small enough neighborhood o/(R"/27rZ") x {0}, k, satisfies 

1^ ■ iy,v) ^ iy + d^lj{r]),ri). 
for some analytic function ip defined in a neighborhood of E C". 

Proof. We may lift k to a mapping between small neighborhoods of R*^ x {0} C C"^^. We 
use the same notations for the lift. By continuity, in small enough neighborhoods we have 
that det(6^) 7^ 0, so by the lemma, locally we may find a analytic generating function ip{y, ^) 
so that the mapping is locally given by 

{y,dy^)^{d^^,0- (11) 

When two neighborhoods overlap, the local generating functions ip, ip must agree modulo 
constants, and in this way we may define a globally defined analytic generating function 
ip{y,C,), so that flTT]) is satisfied everywhere, itself may not give a well-defined function on 
on the base space, but its derivatives are 27rZ-periodic, so flTTl) continues to be valid on the 
torus. 

Still working with the lift, if dy(f = /3(^) is the inverse of the diffeomorphism b, then for 
some function ip, we have 

^{y,0 = m)-y + m 

d^ip = ■ y + dtjj. 

Thus we have 

{y,m))^id^(3-y + di^,0- 
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Because the mapping must be 27rZ-periodic in x and y to descend to the torus, we must 
have that for each ^ C", the assignment Z" 9 n H- (9^/3(^o) defines an automorphism of 
Z". By continuity, d^P G GL(n, Z) must be constant, and since /3(0) = 0, with 9^/3(0) = 1, 
/3 itself must be the identity mapping. 

Therefore, setting ijj = o b{ri), the original mapping is given by 

iy,v) ^ iy + dip, 7]), 

with generating function (p{y, — ^ ' V ~^ i^iO^ ^ is now allowed to be complex- valued. 



Remark 8. Note that Proposition [7] only classifies transformations which preserve indepen- 
dence of the angular variables for a general symbol p{^, e,h) = a-^ + C^(2) (while leaving the 
frequencies a unchanged). That is, if a transformation has this property with respect to any 
such p, then it is easy to see the transformation must satisfy the hypotheses of Proposition 
[71 A priori, some specific p may admit more symmetries than the type described. 

In fact, our main result can be used to prove that this does not happen. Indeed, the normal 
form need not be unique for the reduction process described in the previous section to make 
sense, and then then eigenvalue asymptotics proven in [11] apply. We will see below that 
once we fix a choice of a, by allowing ourselves fiexibility in varying e relative to h, we can 
recover the rest of the asymptotic expansion, which implies uniqueness. Thus we have 

Proposition 9. The quantum Birkhoff normal form of near a Diophantine torus is 
uniquely defined once a choice of the Diophantine frequencies a is fixed. 



4 Eigenvalue asymptotics 

Under the assumptions stated in section [H Theorem 1.1 in [11] implies that if e = 0{h^), 
the eigenvalues of = P -|- ieQ which lie in D have an asymptotic expansion given in terms 
of the QBNF of P^. 

Theorem 10 ([H])- For each N, with p(^) defined above as the "normal form to Nth 
order", let us write 



□ 




(12) 



0(j,k/)<N 




Then for any < 5 < 1, and any fixed integer M, if we require that ^ s < , then 
for each N, as h ^ 0, the quasi- eigenvalues 
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are equal to the eigenvalues of modulo 0{N + 1) in D, in the sense that for sufficiently 
small h, there is a one-to-one partial function from the set of quasi-eigenvalues to spec(Pe), 
equal to 1 -\- 0{N + 1), which is defined whenever a quasi- eigenvalue or the targeted true 
eigenvalue lies in D. 

The expression which appears in place of ^ is the resuh of a Bohr-Sommerfeld type con- 
dition. The constant vector ko contains the Maslov indices and S the actions along a set of 
fundamental cycles of A, for example {xi = 0}, {x2 = 0}. For more details on this point, see 
pT| . as well as section 2 of [TO] . 

5 Main Result 

Our main result is a uniqueness statement roughly stating that if the eigenvalues correspond- 
ing to invariant torus A are the same for operators P + ieQi, P + ieQ2, then they have the 
same QBNF near A. 

Remark 11. As mentioned before, there is a question of whether e may vary independently 
of h. The strongest assumption we consider is that the eigenvalues of two perturbations 
agree for all sufficiently small values of h and for all e such that h'^ < e < (which is 
the range of e- values where the above eigenvalue estimates are valid). In this case we will 
show the QBNF's are necessarily equal. If instead we assume that e = e{h) is some fixed 
function of /i, it is possible that the QBNF's are only equal up to some natural amount of 
symmetry. For example, if e = /i, naturally we cannot extract the coefficients of eH^ and e^h 
from one another. More generally, when e and s is rational we will have the same sort 
of ambiguity. However, as will be explained in the statement of the theorem, up to this sort 
of ambiguity we have uniqueness. 

In the statement of our results we include the possibility that the eigenvalues are only 
known modulo errors which are 0{h^) for some /3. In this case we only recover the terms of 
the QBNF to some order which depends on /3. In general, the presence of noise changes the 
strategy, and as we will see, this indicates that the most robust eigenvalues are those near 
the edge of the window ([3]). 

The main result is the following: 

Theorem 12. Suppose that for v = 1,2, P + ieQy is a pseudodifferential operator satisfying 
the assumptions described in section\^ where P = P'^{x, hD^j.; h) is a fixed, self-adjoint, 
analytic pseudodifferential operator, A is an Hp-invariant Lagrangian torus satisfying the 
Diophantine assumption and Qi = Q"^ , Q2 = Q2 0,1"^ pseudodifferential operators with 
principal symbols qi, q2, respectively, such that (gi) |a = (g'2) |a; and A satisfies the global con- 
dition ^ with respect to each P + ieQy. Let the (e, h)-dependent rectangle D C C be as de- 
scribed in equation with F = (qi) \\ = (^2) |a ■ For u = 1,2, let Pjf^\ denote the QBNF 
of P + i£Qu, which is a formal asymptotic expansion of the form P^T''^ ~ ^jki^ifjkiiO^^^^! 
with P^^li homogeneous of degree j in ^. 
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Suppose that for all sufficiently small h, and e in a range to be specified in cases below, we 
have that 

spec(P + ieQi) = spec(P + zeQg) mod 0{h'^) 

in the window D, in the sense that there is a one-to-one partial function relating the points 
of the two sets which is defined when a point in either set lies in D, and is equal to 1 modulo 
0{hP). 

1. If the above equality holds for all e in the range [h^ , h^], then for each j, k, i such that 

(oo 



{j + k)6 + I < (3, the polynomials P^^ln, u = 1,2 are equal. In particular if (3 = oo 



then all Pjf^^^ agree. 

2. If e = e{h) G [/i^, h^] is some fixed function of h, then we only conclude that for each 
j and for each r G Q such that j5 + r < P, the sums 

,{oo) 



Ep{°°) u - I 2 



(which may be empty, or consist of one or several terms of the expansion) are equal. 
In particular, we have term-by-term equality of the two QBNF's unless for each small 
h, e ^ h'^ for s rational. 

Remark 13. The uniqueness statement above is equivalent to the statement that the QBNF 
can theoretically be recovered by the spectral data, but we will not consider an explicit 
reconstruction process. 

Proof of theorem. We proceed with the proof in two main steps. 
Step 1. Associating eigenvalues to lattice points. 

Given the set of eigenvalues in the window D satisfying the asymptotics f|T2l) . we may not 
know a priori which eigenvalues correspond to which integer lattice points G Z^. However, 
we claim that the real parts of the eigenvalues are enough to determine this correspondence. 
Indeed, the Diophantine condition on a implies that if k, k' G and k ^ k' then 

As /z — 7- 0, by the formula ([5]), the leading term in the real part of the expansion P^°°\h{k — 
^) - £,£,h) is given by 

The values of a, /cq, and S are determined by the (fixed) operator P, so for each k E 7? 
we can calculate the above quantity. (Here we implicitly rely on the choice of action-angle 
variables to make a is well-defined, but it does not matter how we choose.) In the limit 
as /i — )■ 0, this determines the correspondence between eigenvalues in D and lattice points 
keZ\ 
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Step 2. Induction. 



By assumption the QBNF's of P + ieQi and P + ieQ2 agree in all terms which are in- 
dependent of e, since these only depend on the fixed operator P. The strategy is then to 
consider an index (j, k, t) where their QBNF's disagree, which is "minimal" in some sense, 
and to show that the discrepancy must be reflected in the eigenvalues in D. Inductively we 
will then be able to conclude that the two QBNF's agree in at least as many terms which 
do not decay faster than the "noise" term 0{h^). 

Suppose that P[^li 7^ P-fjli some index (j, fc, tj. Then P[^l^ — P2°jM is a homogeneous 
polynomial of degree j which does not vanish identically. By homogeneity, we conclude that 
it is bounded away from in absolute value on some open set f/ C S*^, where denotes the 
unit circle in R^. Then for small enough /z, we can always find ^ = h{k{h) — ^) — ^ such 
that \^\ ~ , with G U. For such ^, we have 

for some C as /;.—)■ 0. Therefore each term Pjfj'keiO^^h^ makes a contribution of the "correct" 
size e^h^^^^ in the spectral asymptotics. 

Say that an index (j, fc, tj such that P[^li 7^ P2^u "minimal" if we can arrange (by 
whatever flexibility we allow ourselves in varying e) that for any other index (j', fc, I') such 
that -Pi°^fc/£/ 7^ P-fji^k't'i '^^ ^s^^Q e^h^^^^ » e^'h^^'^^' as /i — )■ 0. If there exists an index 
(j, /c, €) which is minimal in this sense, and the eigenvalues of P + ieQi and P + ieQ2 agree 
to greater precision than e^h^^~^^, we then have a contradiction. Thus the induction may 
proceed. 

The only issue is that depending on the situation there may not be a single minimal index. 
That indices {j,k,i), {j',k',i') satisfy 

l.,k'j^5f+i' < ^kf^sj+i ^ ce'''h'^'+'' (13) 

for h small enough (for all admissible values of e) defines an equivalence relation on indices, 
and by inductively considering an equivalence class which is minimal in the sense of the 
shared decay rate being the slowest, we may conclude that for each equivalence class / such 
that the common rate of decay is not slower than 0{h^), we have 

Ep(oo) _ p(oo) 

While a priori one might still hope to extricate the individual terms, this does not happen, 
as we will show below. When exactly f|T3|) may hold depends on the cases mentioned in the 
statement of the theorem, so to finish the proof we treat the two cases separately: 

1. If we are allowed the flexibility to vary e to any value between and it is not hard 
to see that we can arrange for e^h^^^^ to decay at a distinct rate for each index (j, /c, t). 
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For example, if we choose e = /i" for some a which is rationally independent of 8 and 1, 
(IT^ never holds, and we get a total ordering of indices. Thus there is always a minimal 
index of discrepancy between p'^^ and ^2°°'*, and the induction shows term- by-term 
that they are equal, for all terms which are not 0{hP\ 

We conclude that P/j^f = ^tfu all (j, k, i) such that 5{j + k) + i < (3. 

2. If e = e{h) is some fixed function of h, it may be that two or more of the decay 
rates satisfy ( fT3l) . However, note that without loss of generality, we can assume that 
whenever (j, k, i) and (j', k', £') satisfy j ^ j', ( |T3i) does not hold. Indeed, when 
choosing the lattice point ^ above, instead of choosing |^| ~ /i'' we could have chosen 
slightly differently, say |^| ~ for S' slightly larger than 6. (fT3|) clearly cannot hold 
for a continuous range of values of 5 if j 7^ j'. 

Thus without loss of generality, the only way f|T3|) can occur is ii < < CK^ 
for some ra, C, i.e. e is comparable to a rational power of /i, such as the case e = h 
mentioned before the statement of the theorem. In this case, we indeed have non- 
uniqueness, as we may add and subtract the same homogeneous polynomial in ^ from 
the coefficients of terms which have the same decay rate, without affecting the eigen- 
values. 

Thus we can only conclude equality of sums Yl Pvjki with the sum taken over each 
equivalence class. As we have seen, equivalence classes do not include indices with 
distinct values of j, so we conclude that for each j G N, r G Q such that 5] + r < (3 
the sums 
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